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We present a low-energy effective field theory to describe the SO(n) bilinear-biquadratic spin 
chain. We start with n = 6 and construct the effective theory by using six Majorana fermions. 
After determining various correlation functions we characterize the phases and establish the relation 
between the effective theories for SO(6) and SO(5). Together with the known results for n = 3 and 
4, a reduction mechanism is proposed to understand the ground state for arbitrary SO(n). Also, 
we provide a generalization of the Lieb-Schultz-Mattis theorem for SO(n). The implications of our 
results for entanglement and correlation functions are discussed. 
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Introduction. - The study of quantum spin chains dates 
back to the early days of quantum mechanics ["]. With 
seemingly simple Hamiltonians, quantum spin chains 
contain very rich physics and thus attract a considerable 
interest. A prominent example is the bilinear-biquadratic 
spin chain with SO(n) symmetry 

h = C o S eY^Y, L f L f+i +s ™ 9 Y,(E L f L f+J 2 - w 

j a<b j a<b 

In the above expression L ab (1 < a < b < n) are the 
generators of SO(n) in the n-dimensional vector repre- 
sentation, with the Casimir operator normalized at ev- 
ery site as J2 a< b(Lj b ) 2 = n— 1. Thus, Eq. (1) represents 
a family of Hamiltonians parametrized by 9 G [— tt, tt] 
and n. As an example, the case of n = 3 corresponds 
to the well-known spin-1 bilinear-biquadratic model, of 
great relevance in the context of Haldane's conjecture [2]. 
Also, for n — 4 the model is equivalent to a symmetri- 
cally coupled spin-orbital chain. Such spin-orbital models 
describe a family of transition metal oxide compounds, 
where orbital degeneracy plays an important role in the 
magnetic properties of the material [3], Quite remark- 
ably, the phase diagram of the model in Eq. (1) for these 
two cases has been established, and exhibits a rich variety 
of phases [4—6]. However, the properties for n > 5 remain 
mostly unclear yet. For instance, it is unknown whether 
the behavior in the large-n limit is somehow similar to 
that of the models with smaller n or not. 

In this context, the main contributions of this Letter 
are (i) to provide a generalization of the Lieb-Schultz- 
Mattis theorem for SO(ra) and (ii) to present an effec- 
tive field theory describing the low-energy physics of the 
model in Eq. (1) for arbitrary n and 9 g [0, #mps], with 
$mps a special point to be discussed later. This is rel- 
evant since it allows us to understand the low-energy 
physics of the system in the considered parameter regime 
in an exact way. To achieve this goal, we start with 
n = 6 and show that at a given phase transition point 
#r, the system is critical and is described by the SO(6)i 
Wess-Zumino-Novikov-Witten model with central charge 
c = 3. This, in turn, allows us to determine various cor- 
relation functions characterizing the phases. Then, after 



establishing a relationship between the effective theories 
for SO (6) and for SO (5) [ ] and gathering some known re- 
sults for n = 3 and 4, we propose a reduction mechanism 
to understand the ground state of the SO(rt) Heisenberg 
chain and suggest an effective field theory for Eq. (1) in 
the considered parameter regime. 

Known results for the general SO(n) case.- Let us now 
comment briefly on what is known about the model in 
Eq. (1) for arbitrary n. In Ref. [8] a phase diagram was 
conjectured based on the existence of some exactly solv- 
able points. For n > 5, this phase diagram exhibits two 
remarkable features between the SO(n) Heisenberg model 
$h = and the integrable SU(n) Uimin-Lai-Sutherland 
model 0tjls = tan -1 [ ]. The first one is a parity 
effect in n; that is, the physical properties of the system 
are sharply different if n is even or odd. This is best seen 
at the special point #mps = tan -1 where the ground 
state of the model is exactly described by a matrix prod- 
uct state (MPS) [8, 10, 11]. This MPS is unique and 
translationally invariant for odd n, whereas it is twofold 
degenerate and breaks translational symmetry for even n. 
The second feature is the location of an integrable model 
at 6* R = tan -1 j^E^p [12], which sits between 6* H and 
$mps for n > 5. This point turns out also to be critical 
for all n, which has important implications. For instance, 
the existence of this point implies that for n > 5 the MPS 
point #mps n0 longer captures the physics of the SO(rt) 
Heisenberg chain. For n = 5 this feature is supported 
by an interesting recent work [7]. Thus, the model in 
Eq. (1) for n > 5 has a quite different phase diagram from 
the n = 3 case, where the Affleck-Kennedy-Lieb-Tasaki 
model at #mps qualitatively describes the properties of 
the spin-1 Heisenberg chain [10]. 

Generalized Lieb-Schultz-Mattis theorem. - From a gen- 
eral perspective, the parity effect has its roots in the dif- 
ference of the SO(n) vector representation for odd and 
even n. More precisely, for even n, there exists an el- 
ement g € SO(n) such that exp(iirg) = —I nX n- The 
presence of this element enables a generalization of the 
Lieb-Schultz-Mattis theorem [1 ]: Assuming that |^) is 
the unique ground state of the model in Eq. (1), the 



2 



e-jtn 




B = -nl2 



FIG. 1: Predicted phase diagram of the S0(6) bilinear- 
biquadratic spin chain in Ref. [ ]. The shaded area is the 
region considered in this Letter. 

"twisted" state |\I> e ) = exp(i-^ Ylf=ij9j)\^) i s n °t oni y 
orthogonal to but also has a vanishing excitation en- 
ergy in the thermodynamic limit N — > oo. This implies 
that the model with even n has cither gaplcss excitations 
or degenerate gapped ground states with broken transla- 
tional symmetry. On the contrary, the model with odd n 
can have a unique gapped ground state. 

Effective theory for n = 6.- Let us consider now the 
model in Eq. (1) for n — 6, whose phase diagram is rep- 
resented in Fig. 1. Since SO(6)~SU(4), this model is 
equivalent to the SU(4) spin chain with self-conjugate 
representation 6 in Ref. [14]. For 9 G [O,0mps], the 
conjectured phase diagram in Ref. [14] is in agreement 
with Ref. [8] and the ground state at #mps is identified as 
an extended valence-bond solid state with broken charge- 
conjugation symmetry. Moreover, an effective theory was 
derived in Ref. [14] by using non-Abelian bosonization 
techniques. Here, though, we use a different approach 
to obtain an effective theory, which makes a clear con- 
nection with the SO(5) effective theory in Ref. [7] and 
enables a possible extension to general n. As expected, 
this effective theory recovers the results from Ref. [ ]. 

Following Ref. [ ], we derive the field theory for the 
SO(6) Heisenberg chain by considering the SU(4) Hub- 
bard model in the strong coupling limit. Let us briefly 
review this approach to make the subsequent discussions 
self-contained and to introduce some necessary notations. 
The SU(4) Hubbard model is written as 

#su(4) - -i^(c] Q c, +1 . Q +H.c.) + [/^(n J -2) 2 , (2) 

3,a j 

where Cj a is the fermion creation operator at site j with 

color index a — 1,...,4 and rij = J2 a c ja c ja i s the 
fermion number operator. Here we consider U > and 
the case of a half-filled fermion energy band (two fermions 
per site). For U >■ t, charge excitation is strongly 
suppressed due to a Mott gap A c ~ U. Thus, two 
fermions frozen at each site constitute six states 4(Cg|0), 
which belong to the vector representation of SO (6). In 
this limit, standard perturbation theory in t/U yields 



an SO(6) Heisenberg model [16]. With this correspon- 
dence in hand, the SO(6) generators L ab are written as 
L ab = lY, a ?cl T Zi3 c f>, where the 4 x 4 SU(4) matri- 
ces T ab are normalized as Tr(T ab T cd ) = 45 ac S bd . Since 
SO(6) is a rank-3 algebra, we choose three diagonal Car- 
tan generators L 12 , L 34 , and L 5e with T 12 = u° <g) a z , 
T 34 = a z ® cr°, and T 56 = a z <g> ct z , where cr° and a z are 
2x2 identity and Pauli matrices, respectively. 

The field theory for the SO(6) Heisenberg chain can 
be derived by applying Abelian bosonization techniques 
to the SU(4) Hubbard model (2) [15]. After lineariz- 
ing the spectra around two Fermi points kp = ±2-7r/ao, 
the fermion operators are decomposed into left-moving 
and right-moving components as Cj a — > y/ao(i>R a e lkFX + 
^L a e~ lkFX ), where a is the lattice spacing. These chi- 
ral fermions are related to boson fields as ipR(L),a = 
(27rao) -1 / 2 Ca exp[±iy / 7r((/) a =p S a )], where Ca are Klein 
factors. The boson fields for U(l) charge and SO(6) 
spin channels are just linear combinations of </> Q , de- 
fined by C) ( ai ) = (0i ± (f>2 + 03 ± 04)/2 and cj) S2 ^ S3 ) = 
(0i ± 02 — 03 T 04)/2 [G, 17], respectively. Similar 
equations hold for their dual fields. The bosonized 
Hamiltonian density % = Hq + Hint contains the 
usual free part H = v -f[K c {d x Q c f + ^(c^) 2 ] + 

2 1 J2l=i[ K s{d x Q Sp ) 2 + Ti- s (dx(t)s p ) 2 } with velocities v c ^ s} 
and Luttinger parameters K c ^ s y The interaction 

part is [15] U ia t = g s Y. P ^ q cos v / 47r0 Sp cos ViTT(f> Sq - 
g sc cos V47r0 c J^l=i cos V47r0 Sp , where g s , g sc > and 
the spin-charge coupling term comes from the 4kp umk- 
lapp scattering presence at half filling. Because of the 
large Mott gap, the charge boson can be safely integrated 
out. By introducing six Majorana fermions to refermion- 
ize the boson fields Sp , for instance (^ + i^ 2 )r(l) ~ 
exp[±iv47i : 0si,ii(i)]) the low-energy effective Hamilto- 
nian density describing the SO(6) spin sector reads [15] 

6 6 

-G s E (3) 

1<I/<K<6 

with m > and G s > 0. Equation(3) predicts a gapped 
ground state with spin-Peierls order for the SO (6) Heisen- 
berg chain [15], which was confirmed numerically [16]. 

Let us now extend the effective theory in Eq. (3) for 
the SO(6) Heisenberg chain in order to incorporate a bi- 
quadratic interaction. In Eq. (3), the Majorana mass 
term is relevant and the fermion interaction term is 
marginally irrelevant for G s > 0. The crucial observa- 
tion here is that, when the Heisenberg chain is perturbed 
by the biquadratic interaction, the only permitted form 
of the effective theory is still Eq. (3) but with modified 
parameters v s , m, and G s as functions of 9. The reason 
for this is that in one spatial dimension only the fermion 
mass term is relevant and the four-fermion interaction 
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term is marginal [18]. Other terms are either irrelevant, 
which play no role in the low-energy limit, or break the 
SO(6) symmetry, which is forbidden since the effective 
theory must preserve the symmetry of the original model. 

More evidence for this effective theory comes from the 
integrable point 6*r = tan -1 | for n — 6. The Bethe 
ansatz solution of this model, as shown by Minahan and 
Zarembo [19], indicates that there are three branches of 
gapless excitations above the SO (6) singlet ground state. 
In the low-energy limit, these excitations have the same 
linear dispersion e(p) ~ v s p. This is in full agreement 
with the effective theory (3) with a vanishing Majorana 
mass, which becomes the SO(6)i Wess-Zumino-Novikov- 
Witten theory (possibly perturbed by marginally irrele- 
vant terms) with central charge c = 6 x | = 3. At this 
critical point, the "speed of light" v s in Eq. (3) is given 
by v s — -|cos0r, which is determined from the Bethe 
ansatz solution. For 9 < Or, the Majorana mass is ex- 
pected to be decreasing when 9 increases from to #r. 
For 9 > #r, it changes sign (m < 0) and an energy gap 
reopens. 

Let us now derive the correlation functions. First, we 
note that the effective theory in Eq. (3) is equivalent to 
six decoupled Ising models and the Majorana mass is 
m ~ (T — T c )/T c with T c the Ising critical temperature 
[20]. For 9 < 9 R , the model in Eq. (1) is in the Ising 
disordered phase with (//„) ^ {v = 1, 2, . . . , 6), while 
for 9 > $r it is in the Ising ordered phase with (cr„) ^ 
(where \i v and a v are the Ising disorder and order oper- 
ators, respectively). In the continuum limit, the SO(6) 
generators are expressed as L ab ~ J^ h + J£ b + (— l) J n ab , 
where the slowly varying SO (6) Kac- Moody currents 
Jr(l) an d the staggered components n ab have critical di- 
mensions 1 and 3/4, respectively. The n ab, s associated 
with the Cartan generators are n 12 ~ 0'i0'2A'3M4A t 5/" t 6! 
n 34 — Mi^203CT4/i5^6, an d n56 ~ ^iM2M3^40"5cr6- Us- 
ing all this information, one can see that the two-point 
correlator (Lj b Lj+ r ) decays exponentially in both Ising 
disordered and ordered phases, while at the critical point 
6>r it decays algebraically as 

(^%)-^ + (-l) r ^ ; ( 4 ) 

where C\ and Ci are nonuniversal constants. Fur- 
thermore, in Ref. [14], a dimerization operator 
v j = {- l )' ] Y,a<b L t L f+i and an operator C 3 = 
{ — ^yj2a<b^j b ^j+i^j+2 characterizing the staggered 
charge-conjugation order were suggested to distinguish 
between the two phases. Remarkably, in the contin- 
uum limit we find that V ~ ^1/^2/^3^4/^5/^6 an d C ~ 
o\<J20-i<J^(j§, satisfying (V,) ^ 0, (Ci) = for 9 < 9 R 
and (T>i) = 0, (Ci) ^ for 9 > #r, which indeed charac- 
terizes the two different phases. At the critical point #r 
these competing orders have the same critical exponents 
and we have (VjV 3+r ) ~ 1/r 3 / 2 and (CjC j+r ) ~ 1/r 3 / 2 . 
Therefore, the critical ground state at #r is an algebraic 



spin liquid unifying the dimerization order, the staggered 
charge-conjugation order, and the Neel order. 

For the solvable point #mpSi it was shown in 
Ref. [8] that the MPS has a hidden antiferro- 
magnetic order quantified by a generalized den Nijs- 
Rommelse string order parameter (SOP) [21] O ab = 
\mi lk _ jHoc (Lfl\^. exp(iTrLf)Ll b ). Because of the 
unbroken SO(6) symmetry, all these SOPs are equal and 
it is sufficient to consider the SOPs for the Cartan gener- 
ators. We find that they are related to the Ising variables 
as 

O 12 ~ (a,) (a 2 > , O 34 ~ (a 3 ) (<r 4 ) , O 56 ~ (a 5 ) (a 6 ) (5) 

These SOPs have nonzero values in the Ising ordered 
phase for 9 > #r and vanish in the disordered phase 
for 9 < #r. Therefore, they are also proper order param- 
eters for the phase with staggered charge-conjugation or- 
der and their utility goes beyond the solvable point #mps ■ 

Reduction from n — 6 to n = 5.- In the following we es- 
tablish a relationship between the effective field theories 
for SO(6) and for SO(5) [ ]. Let us go back to the SU(4) 
Hubbard model in Eq. (2) and interpret the fermion color 
index as spin-| quantum numbers (a = ±|,±i). Then, 
the SO (6) vector representation at each site unifies the 
spin-2 quintet states and the spin-0 singlet state formed 
by two spin- 1 fermions. If an on-site spin-dependent in- 
teraction VJ2 t S 2 is added to the model (2), then the 
singlet and quintet sectors pick up different energies and 
the SU(4) symmetry of the Hamiltonian is broken down 
to SO(5) [22], with an energy difference A s ~ V between 
the two sectors. For V < the quintet states, form- 
ing SO (5) vector representation, have lower energy and 
the effective exchange Hamiltonian in the Mott regime 
is an SO (5) Heisenberg model [22, 23]. In the field the- 
ory treatment [7] , the quintet and singlet degrees of free- 
dom are described by Majorana fermions £ w (y = 1 ~ 5) 
and £ 6 , respectively. In Eq. (3), the energy difference 
A s is formally accounted for by giving £ 6 a large mass 
m s 3> m q > 0, where m q is the mass of the other five 
Majorana fermions. In the low-energy limit, integration 
over £ 6 yields five massive Majorana fermions with a 
renormalized mass m' q > 0. Although the five remain- 
ing Majorana fermions have a smaller mass < m q , 
the SO(5) Heisenberg chain is still in the Ising disordered 
phase, corresponding to dimerized ground states [7]. 

Thanks to this mechanism, we can now study the cor- 
relation functions of the SO (5) model. This is accom- 
plished by simply replacing the Ising variables /x 6 and 
(76 in the expressions for SO(6) by their expectation val- 
ues (ne) 7^ and (er 6 ) = 0, respectively. In the SO (5) 
Ising ordered phase, unlike the SO(6) case, translational 
symmetry is preserved and all two-point correlation func- 
tions decay exponentially. For example, it is easy to show 
that for SO(5) the operator C has a vanishing expectation 
value. However, according to Eq. (5), the nonlocal SOPs 
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are still nonzero and thus are valid order parameters in 
this phase. 

Reduction and effective theory for arbitrary n.- Moti- 
vated by the above discussions, we now propose a reduc- 
tion mechanism to understand the ground state of the 
SO(n) Heisenberg model. As we have seen, the SO(n) 
Heisenberg models for n = 6 and 5 are both described 
by n massive Majorana fermions with mass m > and 
marginally irrelevant terms. Formally speaking, one can 
say that we go from an effective theory for SO(n) to an 
effective theory for SO(n — 1) by giving a large mass to 
the Majorana fermion £ ra and integrating it out. Be- 
cause of the marginally irrelevant couplings in the SO(n) 
effective theory, the elimination of the large-mass Ma- 
jorana fermion induces a negative contribution to the 
mass of the n — 1 remaining Majorana fermions. Usu- 
ally, the exact values of the parameters in these effec- 
tive theories are not known. However, the subtlety here 
is that the SO (4) Heisenberg chain lies exactly at the 
critical point 9 R . Since SO(4)~SU(2)xSU(2), the SO(4) 
Heisenberg chain is equivalent to two decoupled spin-| 
Heisenberg chains, whose effective theory contains four 
massless Majorana fermions [ ]. Thus, if one wishes to 
construct this theory from the effective theory for SO (5), 
the only possibility to obtain the zero mass of the four 
Majorana fermions is the miraculous exact cancellation 
of their original mass [in the SO(5) theory] and the neg- 
ative contribution that comes from integrating over the 
large-mass Majorana fermion £ 5 . We can proceed further 
from this massless SO(4) theory on to an SO(3) theory, 
where three massive Majoranas with m < are obtained 
and which recovers Tsvelik's theory for the spin-1 Heisen- 
berg chain [25]. 

What about the effective field theory for general nl 
By combining the results for n = 3, . . . , 6 and apply- 
ing the reduction mechanism in the reverse direction, 
we can argue that at the integrable point (9r the SO(n) 
spin chain in Eq. (1) is described by an SO(n)! Wess- 
Zumino-Novikov-Witten model with n massless Majo- 
rana fermions perturbed by marginally irrelevant terms. 
This theory has central charge c = n x |, and thus the 
entanglement entropy of a block increases linearly with 
n for the ground state of the system close to and at # R 
[26]. For n > 5, the isolated critical point #r separates 
the region 6 G [0, 6>mps] into and Ising disordered phase 
for 9 < 6*r and an ordered phase for 9 > (?r. By using 
the results in Ref. [27], n zero-energy Majorana modes 
rj a (a = 1 ~ n) localized at the boundary are obtained 
from the SO(n) effective theory for m < in a semi- 
infinite chain and form SO(n) generators T ab = irfn b 
in spinor representation. Note that the SO(n) spinor 
is irreducible for odd n but reducible for even n and 
contains two subspaces. This explains the appearance 
of SO(n) spinor edge states at $mps an d their crucial 
differences for even and odd n [8]. Another interest- 
ing prediction of our theory is the two-point correlator 



(Lj b Lj+ r ) at the critical point 6*r. In the continuum 
limit, the SO(n) generators L ab are a sum of uniform 
SO(n) currents with critical dimension 1 and staggered 
components with critical dimension n/8 (n Ising vari- 
ables). For 3 < n < 7, the contribution from the stag- 
gered components is dominant in the correlator and we 
have {LfLf +r ) ~ (-1)7W 4 . However, for n > 8, the 
contribution from uniform SO(n) currents becomes dom- 



inant and in this case {L ab L ab 
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